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For one flavour, we observe that standard chiral random 
matrix models are schematic variants of the Nambu-Jona- 
Lasinio (NJL) models whether in vacuum or matter. The 
ensuing thermodynamics is that of constituent quarks, with 
mean-field universality in general. For two and three flavours 
non-standard chiral random matrix models with Ua{1) break- 
ing are suggested. For three flavours the transition is remi- 
niscent of the isotropic-nematic transition in liquid crystals. 



1. A large number of finite-temperature and density 
studies in the recent past have focused on the Nambu- 
Jona-Lasinio model (NJL) or extension thereof . These 
models capture the essentials of the spontaneous breaking 
of chiral symmetry in the vacuum, and have been used 
to model the properties of the low-lying pseudoscalar 
nonet. They bear much in common with more mi- 
croscopically motivated thermodynamical descriptions of 
the QCD vacuum such as the instanton model 

In these models, the low-lying spectrum is composed of 
constituent quarks and mesons. In the ground state, the 
quark-antiquark interaction is attractive in the singlet- 
isosinglet providing a simple mechanism for the sponta- 
neous breaking of chiral symmetry, through quark con- 
densation. With increasing density or temperature, the 
constituent quarks in matter overcome the " asymmetry" 
produced by the condensation of the quarks, leading to a 
chirally symmetric phase of screened quarks. The transi- 
tion is mean-field in nature, and mostly entropy driven. 

Since chiral symmetry breaking is encoded in the eigen- 
value distribution at zero virtuality, it is clear that a chi- 
ral phase transition would affect quantitatively this dis- 
tribution. Chiral random matrix models have proven to 
be a very efficient way to get at the eigenvalue distribu- 
tion of QCD inspired Dirac spectra. In this letter, we in- 
vestigate the relationship between chiral random matrix 
models and Nambu-Jona-Lasinio models both in vacuum 
and matter. In section 2, we make this relation quantita- 
tive. In section 3 we discuss the thermodynamics in the 
mean-field approximation. In section 4, we address some 
issues related to the Dirac spectrum in the mean-field 
approximation. In section 5, we discuss new and non- 
standard random matrix models with ?7a(1) breaking as 
inspired by effective models. In section 6, we comment 



on how the issue of confinement may affect the thermo- 
dynamical arguments. Our conclusions are summarized 
in section 7. 

2. To illustrate these points, consider first the case 
of one flavour, with the schematic Lagrangian in four 
Euclidean space as given by 

Ci — +'0^(j7 • d + im + ifijijip 
or equivalently 
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+i^]ii{P + m) Vfl + i^li(P^ + m)tljL + t^PP'' (2) 

in the chiral basis tp = {i^RjipL}- Here P,P^ stand for 
independent auxiliary fields, 5 is a fixed coupling and 
/i is a real chemical potential. Note that the Minkowski 
fields follow from the Euclidean fields through {itp\ ip) 
{ip,^}). Equation (||) is defined on the strip f3 x V3 in 
Euclidean space, with P(t + f3,x) = P{t, x), and ^P{t + 
13, x) = -V'(T,f). 

To establish the connection to chiral random matrix 
models [Q, further simplifications are needed. The anti- 
periodicity of the quark fields yields 
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where ujn — (2n + l)7rT are the Matsubara frequencies 
(T = and x — 1,2,..., N label discrete points 

in space. Space is here a grid of dimension N, where 
each point contains a quark of frequency n. If we fur- 
ther assume that the auxiliary fields P, are constant 
in space and time, the action in (|^) reduces dimension- 
ally to a 0-dimensional one with infinitely many Matsub- 
ara modes. The corresponding partition function can be 
readily found in the form 
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following the rescaling = -/T^V'n and E = V3/2f/^, 
where are now dimensionless Grassmann variables. 
The determinant in (^) is over 2x2 matrices. The cor- 
responding 0-dimcnsional Lagrangian is 

A = + i//)74 + im)q + -^qlqLqnqR (5) 



Consider now the new auxiliary matrix A^'^j with en- 
tries both in ordinary space x, y and frequency space 
n, 771. A is a doubly banded, complex matrix with di- 
mensions {N X N) (g) (oo X oo). In contrast with P, the 
matrix A bosonizes pairs of quarks of opposite chiral- 
ity. For the lowest two Matsubara frequencies it is simply 
{N X N) (g) {2 X 2) matrix . In terms of A, the analogue 
of d) is 



Co = +q^{n + i^)74 + im)q 

+7VETr,,„(AAt) + q^AqL + q^A^ qn 



(6) 



The trace in (|q) is over x and n. The partition function 
associated to p|) is simply 



dAe-^'3STr.,„(AAt)^^^^^^^^^Q (7) 



with the medium Dirac operator in a random back- 
ground, 
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The determinant in (n) is over chirality (2), space (x), 
and frequency space _This is an example of a chiral 
random matrix model 

3. To discuss the thermodynamics of the random 
matrix model, it is best to use (Q) in large N with n — 
iV/Va fixed. The Gibbs free energy associated to (Q) reads 



VL = -N ^w + TlnJ|(l 



(9) 



with Lli'^ — {P + m){P'^ + m). The first term in (^) is 
just the Hartree contribution to the Gibbs free energy 
(Fig. la), while the second term corrects for the double 
counting of the interaction energy (Fig. lb), that is 
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< qlqL >< qnqR > 



(10) 



In this schematic model, the Dirac spectrum is sim- 
plified to two-levels in frequency space (iw) for each 
X = 1,2, N . The quark fields do not carry spin. Recall 
that Np = 1. There is no kinetic energy associated to 
the quarks in either (^) or (^ . These degrees of freedom 
are also manifest in the pressure at high temperature. 



P = = 2nT\ii2-n^PP'' +0{h (11) 
V3 T 

where 2n is the number of quarks and antiquarks, and 
ln2 their respective entropy at T = 00 since their oc- 
cupation number is 1/2. The 1/T term in the entropy 
cancels against the vacuum energy in the Hartree contri- 
bution. The first term in ( pT| ) is the thermal pressure of 
free constituent quarks, while the second term is the left 
out interaction energy. 





(a) (b) 
FIG. 1. (a) is the Hartree contribution to the Gibbs free 
energy, and (b) removes the overcounting. 

In large N, the extremum of the Gibbs free energy (|^) 
yields a gap equation for 

2P*E = l-n-rT (12) 

where n = (e^'^*^^)/^ + 1)^^ for particles and n = 
^g(a),-(-^i)/T_|^-[^-j-i jtqj, antiparticles. The constituent quark 
density is 

PQ = ii'ip^Atp) = -jr^ ^ n{n-n) (13) 
while the quark condensate is 
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.(.^V) = -- = -2nEP. 



(14) 



The relations ([T^ 14) are just schematic versions of the 
usual gap, density and condensate relations in the NJL 
models. 

For zero chemical potential, (|l^) admits solutions for 
4Er* < 1 in the massless case. For zero temperature, 
n = I for fi > Lo and zero otherwise, and n = 0. For 
ji < Lo the pressure of the constituent quarks is constant 
P — n/4E, while that of the free quarks increases linearly 
with /Lt, P = n/i. This behaviour is not compatible with 
confinement at small /i (see section 6). A first order tran- 
sition to the free phase takes place at /x* = u;/2, in which 
chiral symmetry is restored. Clearly the quark number 
density is zero in the constituent quark phase (Fermi level 
in the Dirac gap) and n in the free phase. Since lo plays 
the role of the constituent quark mass, chiral symmetry 
is qualitatively restored for ~ 771^/6, where rriN is the 
nucleon mass. The extra factor of 1/2 follows from the 
peculiar form of the pressure in the absence of kinetic 
energy. 

If we were to retain only few Matsubara modes in the 
random matrix model then the thermodynamical rela- 
tions are altered qualitatively at zero chemical potential 
and quantitatively otherwise. Indeed, at zero chemical 
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potential the Helmholtz free energy with one Matsubara 
mode becomes 
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in comparison with (^) . Setting /? = E ^ to ensure a fi- 
nite T ^ hmit, yields the gap equation u)1 + tt^T^ = 
(jj^/P^E^, for which the pressure is P = — nEP^^. For 
/? = S a single Matsubara mode contributes zero to the 
Hartree part of the free energy (Fig. la). This is rem- 
iniscent of the high temperature cancelation mentioned 
above. We note that due to the choice /? = E the pressure 
grows quadratically in this phase as opposed to linearly 
in (pi]). Since the pressure of a free Matsubara mode 
is Pj;- = 2nln(E7rT)/E, a phase change at ttET* = 1 is 
expected, for which P = P^ = 0. This transition is char- 
acterized by z < '4>'^'4> >— — 2nEP* as an order parameter, 
and was first discussed in |^ . 

In the case where no Matsubara mode is retained and 
/i is finite, the results follow through the substitution 
ttT — > z/i above. In particular, the constituent quark 
mass P* = y/T+'rFJj? /Yi (massless case) grows with ^i. 
Hence i < ip^ip >— — 2nEP* increases in strength, while 
Pq = — 2nE/i decreases. These behaviours are counter- 
intuitive. Indeed, in this phase the pressure P = — nEP^^ 
decreases with increasing fi. Such a system will tend to 
cavitate (collapse) unless a phase transition sets in. This 
is unlikely since the pressure of the free phase is complex 
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The free system is unstable. The present analysis sug- 
gests that the correct thermodynamical behaviour is only 
maintained if all Matsubara modes are retained in the 
random matrix model, with the zero temperature limit 
taken after the large (three) volume limit. 

At zero chemical potential, the random matrix models 
we constructed depend on the number and mixing of the 
various Matsubara modes. At high temperature we ex- 
pect only the lowest modes loq and w_i to survive, mak- 
ing the various matrix models equivalent to each other 
at the critical point. Indeed, at high temperature and 
for T « Tc the gap equation (12) simplifies to 



g{P,+mf + ^I^{P,+m)-h = Oi^) (17) 

with h = m, = {T - Tc)/Tc and g = {12T^)-^ . Eq. 
(|l^) has the generic form of a cubic equation, as expected 
from mean-field treatments of chiral phase transitions . 
Indeed, ( p7| ) is just the gap equation generated from the 
effective potential Oj 
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where in the present case c = 0. This observation im- 
plies that the chiral random matrix model ^ much like 
the NJL model (^) enjoys mean-field critical exponents 
(q, (3, 7, 5, V, 77) = (0, i, 1, 3, i, 0). 

4. To discuss the spectral distribution of the Dirac 
operator in the medium, it is best to use (0-|^), instead 
of (^). For one Matsubara mode, the resolvent of (||) is 



G{z) = — 



Tr 
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(19) 



where the averaging is carried using (Q) in the un- 
quenched approximation. The spectral distribution as- 
sociated to (|l^) follows from its discontinuity through 
the real axis, 7rz^(A) = — ImG(A -I- lO). In the large N 
limit and in the quenched approximation, the resolvent 
( |l9| ) can be readily found to satisfy in general Pastur's 
equation For m = fj, — and at high temperature 

(only one Matsubara mode retained), Pastur's equation 
for G(z) is a cubic equation of the form 



- 2zG^ + (z^ - TT^T^ + 1)G 







(20) 



This equation was first discussed in |^ in a general con- 
text, and used by others [p|,|6| JlO|JTl[ | . The cubic equation 
( pO| ) is a direct manifestation of (|18|) and hence mean-field 
universality H. The spectral distribution at the critical 



point following from (^0|) is consistent with a numerical 
analysis |^ using large matrices. 

For T — m = and finite /i the situation is more 
subtle. Naively, the resolvent in this case satisfies ( |20| ) 
with the substitution ttT i^i. The discontinuities of 
G{z) are valued in the z-plane, with end-points given by 
the zeroes of the discriminant 

V^"* -I- z^iSfi'^ - 20fi^ -1) + 4(^ 1)3 = (21) 

For /X = there are two real roots z = ±2. The result is a 
cut along the real axis between —2 and +2. The discon- 
tinuity along the cut is Wigner's sumrise distribution for 
the spectral density. With increasing /i, the two roots ap- 
proach each other and tend to ±-y/27/8, followed by the 
emergence of two new roots from real infinity. These new 
roots are branch points of cuts that are rejected to infin- 
ity (unbound spectrum). For /ij — 1/8, the branchpoints 
coalesce pair- wise on the real axis, and for /i^ > 1/8 they 
move off the real axis resulting into four endpoints that 
are symmetric about the real axis. The present large N 
analysis suggests that the spectral density jumps to zero 
at /X* = l/\/8E |. Since the pressure ([l^) of the free sys- 
tem is complex, this jump is simply spurious. Indeed, 



^The other choice /3 = 2E amounts for a global shift in the 
pressure, and would not afltect quantitatively the discussion. 



The value of ^, ~ VO.278 quoted in ||6| uses the real part 
of a complex pressure. 
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direct numerical analyses fgjl^l have revealed that the 
spectral distribution at finite /i is complex valued what- 
ever /i > 0. Technically, it is sufficient to note that the 
complex character of the Dirac operator at finite /i causes 
large fluctuations in the eigenvalue distributions, thereby 
upsetting the mean-field approximation for the resolvent 
in certain parts of the z-plane Jl^ , p^ . 

These effects are symptomatic of the fact that the bulk 
pressure shows signs of cavitation (decreasing P) and un- 
stability (complex Pp)- A simple fix would be to include 
all Matsubara modes as we indicated above. Since chiral 
symmetry is only restored for ~ mAr/6 in this case, 
and the thermodynamical system is stable, it would be 
interesting to investigate the properties of the resolvent 
and the Dirac spectrum. 

Quenched lattice simulations ||l^ suggest that the 
transition sets in at /i > 0. This may mean at least 
two thinghs : either that the lattice simulations do not 
generate the proper thermodynamical ensemble, or that 
quenched simulations do not follow from the conventional 
spectrum. The former can be checked by simulating the 
free massive quark ensemble at finite temperature and 
chemical potential, and probing the various limits. The 
latter has been suggested recently in where it was 
argued that quenched lattice simulations follow from a 
new spectrum composed of quarks q and their conjugates 
Q (copies of opposite baryon number). Specifically, the 
analog of (||) is now 



Co 



-q^jiifiq - (5^741^(3 

-^{qliLqRqR + Q[QlQrQr) 



-^{qlQLQRQR + Q\qLq%QR) 



(22) 



As a result, the pressure of the free phase is now real 
= 4nln(E/i)/E and monotonously increasing with /i. 
This pressure is physically meaningful only for S/z > I, 
since a free system with negative pressure will cavitate. 
Besides the chirally broken phase with i < q^q >^ 
and an unphysical pressure (monotonously decreasing 
and negative), (^2|) allows for a chirally symmetric phase 
with a mixed condensate < Q^q > and a physical pres- 
sure. In the large N (Hartree) approximation, the mixed 
condensate satisfies the gap equation 
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(23) 



where r's are Pauli matrices active on the doublet ^ = 
{q,Q)- Hence pi + fi^ = The appearance of 

the mixed condensate, is tantamount to the spontaneous 
breaking of a vector symmetry, that is < ^Vi^ >^ 



Q. The effect of the conjugate quarks is to change to 
— /z^ in the gap equation (thanks to ri) thereby fixing the 
pressure at low ^, 



(24) 



as expected from Fig. lb, since Fig. la contributes to 
zero due to (|2|) for /3 = T,. The quark number density 
in this case is i < C^74T3^ > /V3 = 4nE/z. This be- 
haviour is not compatible with confinement at small fi 
(see section 6). For S/i, — 1, the quark number den- 
sity and the pressure are those of the free constituents, 
so a phase admixture is expected. Using mean-field ar- 
guments, the resolvent associated to (^2|) was discussed 
in j^] in the quenched approximation, and overall con- 
sistency with the numerically generated Dirac spectrum 
was found. 

5. For two and more flavours, we can use the above 
analogy with the NJL model to construct non-standard 
chiral random matrix models with explicit ?/a(1) break- 
ing. Indeed, (|^) suggests the minimal generalization to 
two flavours 



z[T,^l] 



dPdP^ g-W/3E(dct^P-l-dot^pt) 



+ 00 y . , 

n det^,/3(*(™ + ^) 



n— — 00 



(25) 



Here P and are complex 2x2 valued auxiliary fields, 
that transform respectively as VlPV^ and VrP^vI un- 
der U{2)lxU{2)r. In the massless case, (25) is SU{2)lX 
SU{2)fi symmetric. We note that the P-integration can 
be traded to an integration over Grassmann variables, 
leading to a standard two-flavour NJL model with a de- 
terminantal interaction. The QCD motivation for this 
model can be found in |^J^ and references therein. A 
rerun of the above arguments allows a rewriting of ( [25| ) 
into 
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(26) 



where r° = (1, it) and a = 0, 1, 2, 3. For each frequency 
n, the entries in (|2g) are (2, 2) ® {N, N) valued, with m, 

deterministic and Ar^l random and hermitean. 

In the large N limit, the parity even and isospin sym- 
metric saddle point associated to (|2^) hence ( p6| ) belongs 
to the universality class described by the effective poten- 
tial {Np = 2) 



In sau ge t heories this is not in conflict with the Vafa-Witten 
theorem |15| since /x is nonzero. 
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(27) 



is only valid for tempera- 
= 1 — 4ET and for a weak 



with $ - diag Pr.l- Eq.(^ 
tures T r^Tc^ 1 1 AT., with c 
external field h = —to. The quadratic term in (27), orig- 
inates from the 2x2 determinants in ( p5| ) which are seen 
to break explicitly the Ua{^) symmetry. 

For three flavours we can just pursue the analogy with 
one flavour (Q) and two flavours (|2^) and write 



dP ^ptg-A'/3STr/(PP^")+0A^/3(dct_fP+dct/P^") 



+ OC 



n— — oo 



n det^,/3(^("^+^) -^rpn) 



(28) 



where P is a complex 3x3 matrix. In the large N limit, 
the P-integration can be again traded for a Grassmann 
integration, leading to a three-flavour NJL model with 
a determinantal interaction. ( |2^ ) is SU{2>)l x SU{3)ii 
symmetric in the massless case. The QCD motivation for 
( p8| ) can be found in and references therein. In these 
models, we note that the additional parameter 9 > 
plays the role of the instanton density. 

If we were to expand the [/^(l) breaking determinan- 
tal part of the exponential in (|28|), rewrite each deter- 
minant using new auxiliary Grassmann variables (X°, 
i = l,..,N+, and Yj, j = 1,..,7V_, a,b are flavor in- 
dices), and bosonize pairs of opposite chirality as ex- 
plained above (using A's instead of P's) then we can 
rewrite (ESf) in the form 
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N+\N_l \{il3)^ J 



with the (rectangular) random matrix 
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At 
Tl 
\Tr a 
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0^ a 
t / 



(30) 



and the sparse and deterministic (square) matrix 



D 



im r2 + i/i \ 

O + ifj, im 



/ 



(31) 



The block matrices above correspond to the following 
bosonized pairs (the Matsubara indices connected to g's 
are suppressed) 



xy 



1 



3 



(32) 



Each of the four rows and columns in the above matrices 
( ^0[j3l| ) correspond to a vector of length : 1 : 2,® N ® qg, 
2 : ?,®N®oo, 3 : 3®iV+, 4 : 3® A^_, where oo stands for 
the number of Matsubara modes. We observe that ( p9| - 
|3^ with the substitution 3^2, hold for two flavours 
as well, illustrating the model character of the present 
discussion. 

The generic behaviour of these models occur at the 
critical points. For instance, (|2^) undergoes a first or- 
der transition. At the critical point {SU(?>)v symmet- 
ric phase) the universal potential for the parity even 
saddle point with $ ~ diagP follows from (p^, with 
iVj. = 3, /i = -TO, ^2 = 427- - 1, c = A9T > and 
.g = 1/I2r2. In the parameter range {ATT - 1) > and 
{A89^T^-ATT+ 1) > 0, (|l|) allows for a first order tran- 
sition (weak external field) with a critical temperature Tc 
fixed by the conditions £o{Tc,^) = £o{Tc,0) and simi- 
larly for its derivative C'^iTc, $) = C'oiTc, 0). As T ^ 
from below, <i> = 2c/3f; with = 2c^ /%g > 0. For strong 
external fields, the transition turns second order. 

The first order transition described by (|l^) is reminis- 
cent of the one encountered in the Maier-Saupe model 
for the transition from an isotropic molecular phase to 
a nematic phase in liquid crystals [ p^ , pT[ . Indeed, if we 
were to denote by Q the order parameter (mean orien- 
tation of the molecules) in liquid crystals then in the 
mean-field approximation the free energy per molecule is 
easily found to be 



-Tln47r+ (1 -0.4T)Q2 
8 23 4 



105T2 



i75r3 ■ 



(33) 



in analogy with (|18|). The nematic phase corresponds 
to molecules lined up along a preferential direction, but 
with no long-range order. This phase is different from the 
smectic phase which is characterized by various degrees 
of translational ordering. We note that in the presence of 
an external magnetic field (above quark masses) the first 
order transition becomes second order through a tricrit- 
ical point. 

As a final point in this section, we note that in the 
mean- field approximation, the equation of state for 
T ~ allows for a simple understanding of the distribu- 
tion of the quark eigenvalues i'(A) near zero virtuality A ^ 
through the identification 7ri/(A) = Re$(/i ~ iX + Q), in 
the quenched approximation. For two flavours the tran- 
sition is second order with /3 = 1/2, hence = h^^^, 
so that ^(A) ~ A^/3 Por three flavours the transition is 
first order for weak field with i/(0) ~ 2c/ Sir g for T < Tc, 
and i^{0) - for T > T^. 



*From universality and Widom scaling, a similar argument 
would in general suggest ~ A'^/C+t) for the distribution 
of eigenvalues, still in the quenched approximation. Such be- 
haviour, however, calls for a non-local equation of state as 
opposed to a local one in the Landau-Ginzburg analysis. 
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6. The models we have discussed here are schematic 
version of the NJL model and its relatives. For one 
flavour, the spectrum is that of constant a and vr (fluc- 
tuations in P) as well as constituent quarks with masses 
= (P* ~\-Tfi). In conflning theories, constituent quarks 
are barred from the low temperature part of the spec- 
trum. Indeed, in two-dimensional QCD it can be shown 
that the quarks in light-cone gauge satisfy the following 
dispersion relations 



E 



sgn(£' — fc) 



(34) 



where g1 is the fixed QCD coupling in large and A ^ 
is an infrared regulator. For spatially constant modes, 
( p^ ) reduces to -E = ^gH^f^-nX. The spectrum collapses 
to two levels that are rejected to infinity when the in- 
frared regulator is removed. Thermodynamics with the 
spectrum ( p4[ ) yields heavy Boltzmann penalty factors 
g-i/TA ^-p.j\ ^ rpj_^g contribution of the constant con- 
stituent quark modes drop from the thermodynamics at 
low temperature or chemical potential, for A ^ 0. The 
mass gap becomes infinite, and Fig. 1 would be identi- 
cally zero ||2^. As a result, the contribution to the pres- 
sure at low temperature or chemical potential is hadronic. 
So the presently discussed models do not refiect properly 
on conflning theories. They are, however, interesting to 
use around the critical point Tc, since the latter seems to 
be mainly driven by universality. In this regime, there 
may exist a window of coexistence between just about 
to be freed massless quarks and long- wavelength mesonic 
correlations, depending on the nature and character of 
the transition. 

7. We have shown how standard chiral random ma- 
trix models can be understood in terms of two-level NJL 
models, whether in vacuum or matter. The thermody- 
namics of these models is that of constituent quarks, with 
mean-field universality at zero chemical potential. Our 
analysis indicates that the large volume limit does not 
necessarily commute with the zero temperature limit in 
the presence of a finite chemical potential. This point 
may be relevant for lattice simulations of QCD (quenched 
and unquenched) at finite chemical potential llj] • In this 
respect, it would be interesting to reassess the Dirac spec- 
trum at finite chemical potential in the presence of all 
Matsubara modes. 

We have suggested non-standard chiral random matrix 
models for two and three flavours, with explicit Ua^X) 
breaking as inspired by a number of effective models to 
the QCD vacuum. In large N the character of the tran- 
sition follows the general lore of universality |^]. In par- 



^Zero modes are of course an exception. However, in matter 
they do not sense the effect of temperature or density, except 
through a reorganization of the topological configurations. 



ticular for three flavours, and for a speciflc choice of the 
external parameters, the transition is first order and anal- 
ogous to the isotropic-nematic transition encountered in 
liquid crystals. 
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